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' I i The hypothesis of an 'invisible' axion was made by Misha Shifman and 

P i' others, approximately thirty years ago. It has turned out to be an unusually 

^ ' fruitful idea, crossing boundaries between particle physics, astrophysics and 

, cosmology. An axion with mass of order 10~^ eV (with large uncertainties) 

is one of the leading candidates for the dark matter of the universe. It was 
found recently that dark matter axions thermalize and form a Bose-Einstein 
condensate (BEC). Because they form a BEC, axions differ from ordinary cold 
dark matter (CDM) in the non-linear regime of structure formation and upon 
entering the horizon. Axion BEC provides a mechanism for the production 
of net overall rotation in dark matter halos, and for the alignment of cosmic 



microwave anisotropy multipoles. Because there is evidence for these phenom- 
^•-^ , ena, unexplained with ordinary CDM, an argument can be made that the dark 

■ matter is axions. 

o 

O 
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1. Introduction 



It is a great pleasure and honor to be part of Misha Shifman's 60th birth- 
' day celebration. Among Misha's many outstanding contributions to particle 

physics is his well-known proposal, in collaboration with Arkady Vainshtein 
and Valentine Zakharov, that the axion may be very light and very weakly 
coupled [1-3]. Here is the abstract of their paper: 

P- and T-invariance violation in quantum chromodynamics due to the 
so-called 9-term is discussed. It is shown that irrespectively of how the con- 
finement works there emerge observable P- and T-odd effects. The proof is 
based on the assumption that QCD resolves the U(l) problem, i.e., the mass 
of the singlet pseudoscalar meson does not vanish in the chiral limit. We 
suggest a modification of the axion scheme which restores the natural P and 
T invariance of the theory and cannot be ruled out experimentally. 
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The ^-term mentioned by Shifman et al. is 

^« = + ^G'^'^^"''' (1) 

where G^^^, are the QCD field strengths, g is the QCD couphng constant 
and ^ is a parameter. A 0-term is generally present in the action density 
of the Standard Model of elementary particles [4]. Its existence raises a 
puzzle, called the Strong CP Problem. As Shifman et al. explain in their 
paper, the physics of QCD necessarily depends on the value of 9, if none of 
the quark masses vanish, because otherwise QCD wouldn't solve the Ua(1) 
problem (of explaining why the mass of the singlet pseudoscalar meson does 
not vanish in the chiral limit) and hence couldn't be the correct theory of 
strong interactions. This is an important point. If it were possible for QCD 
to be independent of 9, the Strong CP Problem wouldn't be so urgent. 

One can show that QCD physics depends on the value of 9 only through 
the combination 9 = 9 — arg det niq where ruq is the quark mass matrix. If 
9^0 the strong interactions violate P and CP. Such P and CP violation 
is incompatible with the experimental upper bound on the neutron electic 
dipole moment [5] unless |^| < 10"^". In the Standard Model, P and CP 
violation are introduced by letting the elements of the quark mass matrix 
TUq be arbitrary complex numbers [6]. In that case, 9 is of order one. The 
Strong CP Problem is the problem of explaining why |^| < 10^^*^ instead. 

The Strong CP Problem is solved if the term (1) in the Standard Model 
action density is replaced by 

^axion = - ^d^^d^^^ + ^ ^ ^i^-^^""'' 

where (p{x) is a new scalar field, and / is a constant with dimension of 
energy. In the modified theory, 9 = — arg det niq depends on the 
expectation value of <fi{x). This field settles to a value that minimizes the 
effective potential. The Strong CP Problem is solved because the minimum 
of the QCD effective potential V{9) occurs at ^ = [7]. The <fiG ■ G in- 
teraction in Eq. (2) is not renormalizable. However, there is a recipe for 
constructing renormalizable theories whose low energy effective action den- 
sity is of the form of Eq. (2): construct the theory in such a way that it 
has a Upq(1) symmetry which is a global symmetry of the classical action 
density, is broken by the color anomaly, and is spontaneously broken. Such 
a symmetry is called Peccei-Quinn symmetry after its inventors [8]. Wein- 
berg and Wilczek [9] pointed out that a theory with Upq(1) symmetry has 
a light pseudo-scalar particle, called the axion. The axion field is fix), f is 
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of order the expectation value that breaks Upq(1), and is caUed the "axion 
decay constant" . 

The axion mass is given in terms of / by [9] 



All axion couplings are inversely proportional to /. The axion coupling to 
two photons is: 



constant, and is a model-dependent coefficient of order one. It had first 
been thought that / is of order the elcctrowoak scale, in which case the 
axion couplings have strength typical of neutrinos and the axion mass is 
relatively large, in the 10 keV to 10 MeV range. Such axions were quickly 
ruled out by particle physics (beam dumps and rare decays) and nuclear 
physics experiments. But Shifman and others showed [1 3] that / can be 
made arbitrarily large, and hence the axion can be made arbitrarily light 
and weakly coupled. Such an axion is unconstrained by the aforementioned 
experiments, and was dubbed 'invisible'. May an invisible axion really exist? 

The axion has been searched for in many places, and has not been 
found [10]. Axion masses larger than about 50 keV are ruled out by the 
aforementioned particle and nuclear physics experiments. The next range 
of axion masses, in decreasing order, is ruled out by stellar evolution argu- 
ments. The longevity of red giants rules out 200 keV > m > 0.5 eV [11,12] 
in case the axion has negligible coupling to the electron (such an axion is 
called 'hadronic'), and 200 keV > m > 10~^ eV [13] in case the axion has 
a sizable coupling to electrons. The duration of the neutrino pulse from 
supernova 1987a rules out 2 eV > to > 3 • 10~^ eV [14]. Finally, there is a 
lower limit, m ^ 10~^ eV, from cosmology which is discussed in the next 
section. This leaves open an "axion window" : 3 • 10~^ > 10~^ eV. The 
lower edge of this window (10~^ eV) is much softer than its upper edge. 

2. Axion production in the early universe 

There are two populations of axions produced in the early universe, which 
we may call 'hot' (or thermal) and cold. Hot axions are produced in thermal 

processes such asg + y— >(7 + a where g is a quark and g a gluon, or 
TT + TT — * TT + a where tt is a pion [15-17]. The number density of thermal 



m ~ 6 eV 



10^ GeV 
/ 



(3) 
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axions today (time to) is 



where A/d is the effective number of thermal degrees of freedom at the time 
axions decouple from the thermal bath. The Standard Model has Ad = 
106.75. Thermal axions are a form of hot dark matter, similar to neutrinos, 
in the context of large scale structure formation. 

The cold axions are produced when the potential V^(^), and hence the 
axion mass, turns on near the QCD phase transition [18]. The critical time, 
defined by m(ii)ti = 1, is ti ~ 2 • 10"'^ sec (//10-^^GeV)3 . Cold axions are 
the quanta of oscillation of the axion field that result from the turn on of 
the axion mass. The average number density of cold axions at time ii is 

n„(ii) ^m{h){^^h)) ~ TT/^i . (6) 

In Eq. (6), we used the fact that the axion field cp{x) is approximately 

homogeneous on the horizon scale ti, because wiggles in (p{x) which entered 
the horizon long before ti have been red-shifted away [19]. We also used the 
fact that the initial departure of tp{x) from the CP conserving minimum is 
of order /. The axions of Eq. (6) are non-relativistic. Assuming that the 
ratio of the axion number density to the entropy density is constant from 
time ti till today, one finds [18] 

"-Kio4v)'(¥)' 

for the ratio of the axion energy density to the critical density for closing 
the universe, h is the present Hubble rate in units of 100 km/s.Mpc. The 
requirement that axions do not overdose the universe implies the constraint 
m > 6 • 10~^ eV. For a more detailed discussion of the production and 
properties of dark matter axions, the reader may wish to consult refs. [17, 
20]. 



3. Dark matter caustics 

It has been established from a variety of observational inputs that approx- 
imately 23% of the energy density of the universe is "cold dark matter" 
(CDM). The CDM particles must be non-baryonic, and cold. "Cold" means 
that their primordial velocity dispersion is small enough that it can be set 
equal to zero for all practical purposes when discussing the formation of 
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large scale structure. The leading candidates for the CDM particles are 
axions, weakly interacting massive particles (WIMPs), e.g. the neutralino 
in supersymmetric extensions of the Standard Model, and sterile neutrinos 
with mass in the keV range. 

A central problem in dark matter studies is the question how CDM is 
distributed in the halos of galaxies, and in particular in the halo of our 
Milky Way galaxy. Indeed, knowledge of this distribution is essential for 
understanding galactic dynamics and for predicting signals in direct and 
indirect searches for dark matter on Earth. 

Galactic halos are thought to be collisionless fluids and must there- 
fore be described in 6-dimensional phase space. A full description gives the 
phase space distribution f{r,v;t) of the dark matter particles in the halo, 
i.e. their velocity (v) distribution at every position f. An important sim- 
plification occurs in the case of cold dark matter because CDM particles 
lie in phase space on a thin 3-dimcnsional hypersurface. This fact implies 
that the velocity distribution is everywhere discrete [21] and that there are 
surfaces in physical space, called caustics, where the density of dark matter 
is very large. 

Calactic halos have two types of caustics, outer and inner. The outer 
caustics are simple fold {A2) catastrophes located on nested topological 
spheres surrounding the galaxy. The catastrophe structure of the inner 
caustics depends on the angular momentum distribution of the infalling 
dark matter particles. If that angular momentum distribution is dominated 
by net overall rotation, implying V x w 7^ 0, the inner caustics are a set of 
'tricusp rings'. A tricusp ring is a closed tube whose cross section is a sec- 
tion of the elliptic umbilic (-D-4) catastrophe [22 24]. The rings are located 
in the plane of the galaxy. In the self-similar infall model [25], generalized 
to include the effect of angular momentum [26,27], the caustic ring radii a„ 
(n = 1, 2, 3..) are predicted to obey the law [22,27] 



where iij-ot is the rotation velocity of the galaxy and jmax is a parame- 
ter characterizing the amount of angular momentum that the dark matter 
particles carry. 

Observational evidence for caustic rings at the radii a„ predicted by Eq. 
(8) was found in the rotation curves of external galaxies [28] , the rotation 
curve of our own galaxy [29], and an IRAS map of the Galactic plane in the 
direction of the nearest caustic ring [n = 5) [29]. A summary of the evi- 
dence can be found in ref. [27] . The evidence implies that the distribution of 




(8) 
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imax values over nearby spiral galaxies, including the Milky Way, is peaked 
at jmax — 0.18. It also implies that we on Earth are close to a cusp in the 
nearest caustic ring of dark matter. As a result, the dark matter velocity 
distribution on Earth is dominated by a single flow, of known velocity vec- 
tor. That single flow, called the "Big Flow" has density of order 1 GeV/cc, 
which is two or three times larger than the commonly cited estimates of 
the total local dark matter density. 

Finally, the evidence for caustic rings of dark matter halos implies that 
the dark matter particles fall in with net overall rotation, and hence that 
their velocity field has non-zero curl: V x tf 7^ 0. If their velocity field 
were irrotational (V x v = 0), the inner caustics would have a tent-like 
structure [24] which is quite distinct from that of the caustic rings for 
which evidence was found. 

This raises a puzzle. Indeed if the dark matter is cold and collisionless, 
as is the case for weakly interacting massive particles (WIMPs) and was 
thought to be the case for axions, the velocity field remains irrotational at 
all times because it is the outcome of gravitational forces proportional to the 
gradient of the Newtonian potential [24] . (General relativistic effects allow 
the creation of rotational velocity fields but are subdominant because the 
velocities involved are much less than the speed of light.) Thus, if the dark 
matter is cold and collisionless, one expects the inner caustics of galactic 
halos to be the tent-like structures of the W xv — case, instead of the rings 
for which evidence was found. This puzzle has bothered me for a number of 
years, but Qiaoli Yang and I may now have found a solution to it [30]. As 
explained below, cold dark matter axions form a Bose-Einstein condensate. 
As a result, their properties differ from those of ordinary CDM. 

4. Bose-Einstein condensation of dcirk matter axions 

a 

The number density of cold axions implied by Eq. (6) is 
, , 4-1047 / / \ t /a(ti)\3 

where a(t) is the cosmological scale factor. Because the axion momenta are 
of order at time ti and vary with time as a{t)~^, the velocity dispersion 



^'AU the material in this section is taken from ref. [30]. 
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Sv{t) ^ J- ^ (10) 
^ ^ mh a{t) ^ ' 

i/each axion remains in whatever state it is in, i.e. if axion interactions are 
negligible. Let us refer to this case as the limit of decoupled cold axions. If 
decoupled, the average state occupation number of cold axions is 

Clearly, the effective temperature of cold axions is much smaller than the 
critical temperature 

T = " ~ 300 GeV ( ^ \ " ^ (12) 

for BEC. Axion number violating processes, such as their decay to two 

photons, occur only on time scales vastly longer than the age of the universe. 
The only condition for axion BEC that is not clearly satisfied is thermal 
equilibrium. 

Axions are in thermal equilibrium if their relaxation rate T is large 
compared to the Hubble expansion rate H(t) = ^. At low phase space 
densities, the relaxation rate is of order the particle interaction rate Tg = 
na6v where a is the scattering cross-section. Axions have self interactions 
described by the action density /!scif = +\^'P'^ where A ~ 0.35(y)2. The 
cross-section iov ip + (f ^ (f + (f scattering due to axion self interaction is 
in vacuum 

If one substitutes (Tq for a, F., is found much smaller than the Hubble rate, 
by many orders of magnitude. However, in the cold axion fluid background, 
the scattering rate is enhanced by the average quantum state occupation 
number of both final state axions, a ^ (T()J\f'^, because energy conservation 
forces the final state axions to be in highly occupied states if the initial 
axions are in highly occupied states. In that case, the relaxation rate is 
multiplied by one factor of J\f [31] 

T'^naodvAf . (14) 

Combining Eqs. (9-11,13), one finds T{ti)/H{ti) ~ 0(1), suggesting that 
cold axions thermalize at time ti through their self interactions, but only 
barely so. 
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It may seem surprising that the huge and tiny factors on the RHS of 
Eq. (14) cancel each other. In fact the canceUation is not an accident. Con- 
sider a generic axion-like particle (ALP) whose mass m and decay constant 

2 

/ are unrelated to each other. Its self interaction coupling strength A ~ y^. 
Cold ALPs appear at a time ti ^ ^ with number density n{ti) ^ f^m, 
and velocity dispersion 6v{ti) ~ 1. Substituting these estimates in Eqs. 
(11), (13) and (14), one finds that the thermalization rate is of order the 
Hubble rate at ti, for all / and m. 

A critical aspect of axion BEC phenomenology is whether the BEC con- 
tinues to thermalize after it has formed. Axion BEC means that (almost) all 
axions go to one state. However, only if the BEC continually rethermalizes 
does the axion state track the lowest energy state. 

The particle kinetic equations that yield Eq. (14) are valid only when 
the energy dispersion ^m{Sv)'^ is larger than the thermalization rate [31]. 
After ti this condition is no longer satisfied. One enters then a regime where 
the relaxation rate due to self interactions is of order 

Fa - A n m'^ . (15) 

T\{t)/ H{t) is of order one at time ti but decreases as t a{t)~^ afterwards. 
Hence, self interactions are insufficient to cause axion BEC to rethermalize 
after t\ even if they cause axion BEC at ti. However gravitational inter- 
actions, which are long range, come in to play. The relaxation rate due to 
gravitational interactions is of order 

Pg ~ G n TO^ £2 (^g^ 

where ^ ~ {m5v)~^ is the correlation length. Tg{t)/H{t) is of order 
4 • 10-^(//10i2 GeV)t at time t\ but grows as ta-'^{t) oc a{t). Thus gravi- 
tational interactions cause the axions to thermalize and form a BEC when 
the photon temperature is of order 100 eV (//lO^^ GeV)2. 

The process of axion Bose-Einstein condensation is constrained by 
causality. Thus one expects overlapping condensate patches with typical 
size of order the horizon. As time goes on, say from t to 2t, the axions in 
t-size condensate patches rethermalize into 2i-size patches. The correlation 
length is then of order the horizon at all times, implying iJf ~ instead 
of Eq. (10), and T^/H cx f^a-'^{f) after the BEC has formed. Therefore 
gravitational interactions rethermalize the axion BEC on ever shorter time 
scales compared to the age of the universe. 

The axion field may be expanded in modes labeled a: 

Vix) = ^ [as ^s{x) + a], (17) 

a 
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where the $5(3;) arc the positive frequency c-number solutions of the 
Heisenberg equation of motion for the axion field 

D^D,^{x) = ,9^'19^a, - Vl,d>]ip{x) = m\{x) , (18) 

and the and at are creation and annihilation operators satisfying canoni- 
cal commutation relations. We neglect the self-interaction term which would 
otherwise appear on the RHS of Eq. (18), because it is of order -p'fi where 

p is the axion density, and hence smaller by the factor {^^^fy^ than the 
relevant terms (of order ™ equation. BEC means that all cold 
axions, except for a small fraction, go to a single state which we label (5 = 0. 
The corresponding ^o{x) is the axion wavefunction. In the spatially flat, 
homogeneous and isotropic Robertson- Walker space-time, 

a{t)2 

where A is a constant. The state of the axion field is >= 
(I/a/M) (ao)^|0 > where |0 > is the empty state, defined by as |0 > 
= for all a, and A'' is the number of axions. 

To compare axion BEC with CDM, let us divide the observations into 
three arenas: 1) the behaviour of density perturbations on the scale of the 
horizon, 2) their behaviour during the linear regime of evolution within the 
horizon, and 3) their behaviour during the non-linear regime. CDM provides 
a very successful description in arena 2. However, axion BEC and CDM are 
indistinguishable in arena 2 on all scales of observational interest [30,32]. In 
particular, the equation governing the evolution of axion BEC perturbations 
is 

d^5 + 2Hdt6 - (^AttGpo - S = (20) 

where k is co-moving wavevector. The last term in Eq. (20) is absent for 
CDM. Eq. (20) implies that the axion BEC has Jeans length 

fc-i = (167rGpm^)-^ 

=i.o2.io-cmri^yr^^^^^^V . (^d 



m J \ p 

However, the Jeans length is small compared to the smallest scales (^ 100 
kpc) for which we have observations on the behavior of density perturba- 
tions in the linear regime. 

In the non-linear regime of structure formation (arena 3) and in the 
absence of rethermalization, axion BEC and CDM again differ only on 
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length scales smaller than de Broglie wavelength. This follows from the 
WKB approximation and has also been shown by numerical simulation [33] . 
Since the axion de Broglie wavelength (of order 10 meters in galactic halos) 
is negligbly small compared to all length scales of observational interest, 
we again find that axion BEC and CDM are indistinguishable when there 
is no rethermalization of the BEC. 

However, it was found above that gravitational interactions do rether- 
malizc the axion BEC continually so that the axion state tracks the lowest 
energy state. This is relevant to the angular momentum distribution of 
dark matter axions in galactic halos. The angular momentum of galaxies is 
caused by the gravitational torque of nearby galaxies early on when pro- 
togalaxies are still close to one another [34]. As was mentioned in Section 
III, CDM presents us with a puzzle. The velocity field of ordinary cold 
dark matter, such as WIMPs, remains irrotational whereas the evidence 
for caustic rings of dark matter implies that the dark matter falls in with 
net overall rotation. The puzzle is solved if the dark matter is an axion 
BEC which rethermalizes while tidal torque is applied to it. Indeed, the 
lowest energy state for given total angular momentum is one in which each 
particle carries an equal amount of angular momentum. In that case there 
is net overall rotation. V x w 7^ is accomodated in the BEC through the 
appearance of vortices. The phenomenon is observed in quantum liquids 
and well understood [35]. 

Finally let's consider the behaviour of density perturbations as they 
enter the horizon (arena 1). Here too axion BEC differs from CDM. The 
CDM perturbations evolve linearly at all times. The axion BEC perturba- 
tions do not evolve linearly when they enter the horizon bccaiisc the conden- 
sates which prevailed in neighboring horizon volumes rearrange themselves, 
through their gravitational interactions, into a new condensate for the ex- 
panded horizon volume. This produces local correlations between modes of 
different wavevector since the perturbation of wavcvcctor fc, upon entering 
the horizon, is determined by the perturbations of wavevector say ^A; in its 
neighborhood. We propose this as a mechanism for the alignment of CMBR 
anisotropy multipolcs [36] through the integrated Sachs- Wolfe (ISW) effect. 
Unlike CDM, the ISW effect is large in axion BEC because the Newtonian 
potential changes entirely after entering the horizon in response to the 
rearrangement of the axion BEC. 
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